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Abstract. We give three formulas for meromorphic eigenfunctions (scat- 
tering states) of Sutherland's integrable A'^-body Schrodinger opera- 
tors and their generalizations. The first is an explicit computation of 
the Etingof-Kirillov traces of intertwining operators, the second an in- 
■ tegral representation of hypergeometric type, and the third of Bethe 

ansatz type. The last two formulas are degenerations of elliptic formulas 
obtained previously in connection with the Knizhnik-Zamolodchikov- 
Bernard equation. The Bethe ansatz formulas in the elliptic case are re- 
^ ' viewed and discussed in more detail here: Eigenfunctions are parametrized 

, by a "Hermite-Bethe" variety, a generalization of the spectral variety 

of the Lame operator. We also give the g-deformed version of our first 
formula. In the scalar sIn case, this gives common eigenfunctions of the 
commuting Macdonald-Rujsenaars difference operators. 
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1. Introduction 



Let g be a simple complex Lie algebra, with a non-degenerate ad-invariant 
Q I bilinear form ( , ), and a fixed Cartan subalgebra f) . Let A C f) * be the set 

' of roots of g . For each a € A let G g be a corresponding root vector, 

normalized so that (ea,e_Q,) = 1. 

Suppose that C/ is a highest weight representation of g with finite dimen- 
sional zero- weight space ( the space of vectors in U annihilated by f) ) U[0]. 



\ We consider in this paper the differential operator 

2 

^ = -A + ^ -^/^----e„e_„ (1) 
^sm^(7ra(A)) 

acting on functions on f) with values in U[0]. The Laplacian A is the op- 
erator X^j, 7^ in terms of coordinates = {bu,X) for any orthonormal 
basis bi, . . . ,br of . As remarked in ||3[, if g = sIn, and U is the sym- 
metric power S^^C^ of the defining representation C^, then U[0] is one- 
dimensional and this differential operator reduces to Sutherland's integrable 
A^-body Schrodinger operator [Iq] in one dimension with coupling constant 
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p{p +1). In suitable variables, H is in this case proportional to 

Hs = -T.3Z2+PiP+^)Yl Vi^^ - ^i) (2) 

with V{x) = 1/ sin^(x) (trigonometric model) or V{x) = 1/ sinh^(a;) (hy- 
perbolic model). We refer to this special cases as the scalar cases. 

We consider in this paper eigenfunctions (functions ip with Hip = eip for 
some e G C) of the form e^'^*^^''^V(A) where / is meromorphic and regular 
as a(A) too, Va G A (a more precise definition is given below). In the 
language of A^-body Schrodinger operators these are scattering states for the 
hyperbolic models. It turns out that for generic the space -E'(^) of such 
eigenfunctions is finite dimensional and isomorphic to U[0]. 

We give three formulas for these eigenfunctions. The first one (Theorem 



is an explicit cumputation of an expression of Etingof and Kirillov in 
terms of traces of certain intertwining operators. The reason that these 
combinations of traces can be computed explicitly is that they turn out 
to be the same as for the corresponding Lie algebra without Serre rela- 
tions (Prop. |9.2[ ). This observation reduces the computation of traces to 



a tractable combinatorial problem. We also give the g-deformed version of 



this formula (Theorem 8.1). The second formula, an integral representation 
(Theorem fl.l[ ) , follows from the fact that eigenfunctions can be obtained as 
suitable limits of solutions of the Knizhnik-Zamolodchikov-Bernard equa- 
tion of conformal field theory, for which we gave explicit solutions in Q . The 
third formula (Theorem 5.1) is of the Bethe ansatz (or Hermite) type: one 



has an explicit expression of a function depending on parameters T € C". 
This function is an eigenfunction if the parameter lie on an algebraic vari- 
ety which we call Hermite-Bethe variety. We have a regular map p from 
the Hermite-Bethe variety to f) , sending T to the corresponding value of 
^. Conjecturally, p is a covering map and a basis of i?(^) for generic ^ is 
obtained by taking the eigenfunctions corresponding to the points in the 
fiber p~^{C)- This conjecture is proved in some cases including the scalar 



case (Theorem 5.3). 

The three formulas hold for generic values of the spectral parameter ^. 
In Section ^ we study trigonometric polynomial solutions of the eigenvalue 
problem Hip = eip. These (Weyl antiinvariant) eigenfunctions are related 
to multivariable Jacobi polynomials (or Jack polynomials). Formulas for 
these polynomials are obtained as the spectral parameters tends to an an- 
tidominant integral weight. The construction involves the construction of 
the scattering matrices for the problem, which give a solution of the Yang- 
Baxter equation. It would be interesting to generalize this construction to 
the g-deformed case, which would give formulas for Macdonald polynomials. 
We give such formulas in Section ^ as a conjecture. This conjecture can be 
proved by the same method as in the classical case for sl2 and but the 
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regularity property of Weyl antiinvariant eigenfunction appears to be more 
difficult to prove in the general case. 

The fact that we have three different formulas for the same thing can be 
explained in informal terms as follows. Our second formula is an integral 
depending on a complex parameter k ^ 0. Its form is 

V;(A) = e2^^«(^) / $o(r)^/''tj(r,A) 

The function is a many-valued homolomorphic function on Dn = (C — 
{0, 1})'^ — Uj<j{T I Ti = Tj} and w(T, A) is a rational ?7[0]-valued differential 
n-form on Dn- The cycle 7(k) has coefficients in the local system determined 
by ^o{T)^^'^- It turns out that, for all k, 7(k), ip is an eigenfunction of H 
with the same eigenvalue 47r^(^,^), and that (if is generic) for any fixed 
generic k, all eigenfunctions can be obtained by suitable 7(k). One can thus 
consider suitable families of cycles parametrized by k in the limits /t ^ 0, 
K ^ oo. In the former case the homology reduces to ordinary homology, 
and the answer is given in terms of residues of uj (our first formula) . In the 
latter limit, the integral can be evaluated by the saddle point method and 
the eigenfunctions are given by evaluating uj at the critical points of (our 
third formula). 

This reasoning turns out to be useful to write down the formulas, (and to 
"understand" them) but it is then easier to prove them by other methods. 

In Section ^, we generalize our results on Bethe ansatz eigenfunctions (the 
third formula) to the elliptic case. The Schrodinger operator in this case has 
the Weierstrass p-function as a potential: 

He = —A + p(a(A))eQ,e_Q, + const. (3) 

In the scalar case q = sl]\j, U = 5P^C^, we obtain (|2|) with V = p. 
Bethe ansatz eigenfuctions of Hg were given in ||^]. After reviewing this 
construction, we give a result on completeness of Bethe states (Theorem 
[7. 3D parallel to the trigonometric case, in the scalar case and the case of the 
adjoint representation of sIat: for generic values of the spectral parameter ^ 
there exist dim(C/[0]) solutions of the Bethe ansatz equations corresponding 
to linearly independent quasi-periodic eigenfunctions with multiplier given 

The formulas discussed here in the trigonometric case were obtained by 
considering the trigonometric limit of solutions of the KZB equations for 
an elliptic curve with one marked point. The same construction could be 
done in the case of N marked points. It turns out that in the trigonometric 
limit one of the equations always reduces to an eigenvalue equation for H. 
The new feature is that U is the tensor product of N representations. The 
remaining commuting operators define differential equations in the space of 
eigenfunctions of H with a fixed eigenvalue. 
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In the g-deformed case, we could only generalize the first formula. It 
would be interesting to find the second and third formula in the g-deformed 
case. 

The organization of the paper is as follows: in Section ^ we define and 
describe the space of eigenfunctions we consider in the trigonometric case. 
Three formulas for these eigenfunctions are given in Section |3| in terms of 
coordinates of singular vectors, in Section ^ as integrals, and in Section |5| by 
a Bethe ansatz. In Section |6|, the action of the Weyl group is discussed along 
with the relation to multivariable Jacobi polynomials at special values of the 
spectral parameter. In Section ^ we discuss the Bethe ansatz in the elliptic 
case. Section ^ gives a generalization of the first formula to the (/-deformed 
case. The proof of the first formula is contained in Section p. 

We conclude this introduction by fixing some notations and conventions. 
The Cartan subalgebra f) will be often identified with its dual via ( , ). We 
fix a set of simple roots qi, . . . ,0^ € A, and write Q for the root lattice 
©jZoj. Its positive part ©jNoj, with N = {0, 1, . . . }, will be denoted by 
Q+. We denote by p half the sum of the positive roots. We will use the 
partial ordering /?> /3' iff /?-/?' G on f} and write (3 > (3' \i (3 > f3' but 
(3 / 13'. We set A+ = A n Q+. If /3 € Q, we write \I3\^ = {(3,13). We will 
also use the notation \A\ to denote the cardinality of a set A. The group of 
permutations of n letters is denoted by Sn- 



2. The il^ functions 

We introduce a space of meromorphic functions on \] that is preserved by 
the differential operator H. For /3 G f) *, let 

X/3(A) = e-2-^W, (4) 

and write Xp = Xi \i (3 = Ui. Let A be the algebra of functions on \) 
which can be represented as meromorphic functions of (^1, . . . ,^r) £ 
with poles belonging to set UogAj^a = For instance, if a is a positive 
root, the functions 

■ 2/ (^^^ = ~ n ^'^v cot(7ra(A)) = i |^^° (5) 

sm (7rQ(A)) [I — X^y 1 — X^ 

belong to A. The algebra ^4 is a subalgebra oi A = C[[Xi, . . . , Xr]] 

If ^ € we introduce the A-module A{(^) and the 74-module A(^) of 

functions of the form exp(27ri^)/ where f & A, resp. f £ A. These modules 

are preserved by derivatives with respect to A. Therefore H preserves the 

spaces A{0 U[0], A{0 U[0]. 

Thus any ip € A{(^) ® U[0] has the form 

/3eQ 



with ijjp G U[0] vanishing if /? Q+. 
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Theorem 2.1. For generic ^ € f) *, and any non-zero u G U[^], there exists 
a unique ip = YljS -^13-^^13 ^ ^iO ^ ^^[0]; such that Hip = exp, for some 
€ G C, and such that ipo = u. Moreover, e = {2it)'^ {S^ , , andip € ^(^)(8C/[0]. 

The existence and uniqueness of as a formal power series was proved in 
the scalar case by Heckman and Op dam (see Sect. 3 of [^), who generalized 
a classical construction of Harish-Chandra (see IV. 5). The fact that the 
series converges to a meromorphic function seems to be new. 

We next prove this theorem except for the statement that the formal 
power series ip actually belongs to A{S^) (8) U[0], which will follow from the 
explicit formulas for ip given below. 

Proof: The idea is that the eigenvalue equation Hip = eip is a recursion 
relation for the coefficients t/)^. We use the fact that a and —a give the same 
contribution to the sum in H, to apply (I). With the formula AXp_^{X) = 
(27ri)^(/3 — — C)Xp-^{X), we see that Htp = eip is equivalent to the 
recursion relation 

[(/3 -CP-O- (2vr)-2e]V;^ = ^ eae-ai'p-ja (6) 

for the coefficients V/S £ ^[0]- The sum on the right-hand side has only 
finitely many non-zero terms. The initial condition for this recursion is 
ipo = u 0. For /3 = 0, the equation reads then {CO ~ (2vr)~^e = 0. Thus 
there is a solution only if e = (27r)^(^, 0- For generic ^, the coefficient of ipp 
does not vanish if /3 / 0, so d6|) gives ipp in terms of ippi with /3' < /3. We 
conclude that, for generic ^, (^ has a solution if and only if e = (27r)^(^,^), 
and this solution is unique. □ 

Definition: We denote by E{0 the vector space of functions ip G A(^)(g)f7[0] 
such that Hip = {2-Kf{^,0^. 

Remarks. 1. It looks as if the definition of E{£P) should depend on the 
choice of simple roots, but this is not so. From the explicit expressions 
given below it follows that, for any set R of simple roots, the functions in 
E{^) can be analytically continued to functions in A/j(^)(gi?7[0], where Aii{SP) 
is the space A{0) defined using R. See Section ^ 

2. if is part of a commutative r-dimensional algebra D of Weyl- invariant 
differential operators whose symbols are (Weyl-invariant) polynomials on 
\)* . These operators have coefficients which are polynomials in cot(7ra(A)), 
a G A_(_. It follows that A{£P) is preserved by see (^), and that the 
functions ip of the theorem are common eigenfunctions of all operators in D. 
Indeed \i L G D with symbol P, then L preserves E{SP) since it commutes 
with H. If G E{0 has leading term e'^'^^^^^^u, u G C/[0], then Lip has 
leading term P(27ri^)e2'^*?(^)u. Thus Lip = P{2mCjip. 

3. If g is a general Kac-Moody Lie algebra and C/ is a g -module with finite 
dimensional zero- weight space f7[0], H is still a well-defined endomorphism 
of A{£P) ® ?7[0]: the (possibly infinite) sum over A_|_ gives only finitely many 
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contributions in each fixed degree, as sin^^(7ra(A)) = 0{Xa), see @. The 
above theorem holds (except for the statement that V ^ ^(C) ® U[0]) with 
the same proof. 

Etingof and Kirillov Q gave a representation theoretic construction of 
eigenf unctions: given ^ generic and u G U[0], let M^_p be the Verma module 
with highest weight C ~ P- Then there exists a unique homomorphism G 
Homg (M^_p, Mg_p (g) ?7), such that the generating vector v^_p € M^_p is 
mapped to f^-p u + • • • , up to terms whose first factor is of lower weight. 
Then ip € ® U[0] is the ratio of formal power series 

tr/ift exp(27riA) 

M>^) = l'^' ■ (7) 

trM_p exp(27r?A) 

The traces are formal power series whose coefficients are traces over the 
finite dimensional weight spaces of the Verma modules. By definition, the 
trace of a map M ^ M ®U is the canonical map (for finite dimensional M) 

tTM ■■ Homc(M, M U) ^ M* M U ^ U. 

The numerator in (|^ belongs to — p) (>S> U[0], and the denominator to 
A{—p) with leading coefficient 1, so the ratio is a well-defined element of 



U[0]. Combining this result with Theorem 2T, we get 



Proposition 2.2. For any generic ^ G f) *, the map u ^ tpu defined by (Q) 
is an isomorphism from U[0] to i?(^). 

Remark. Etingof and Styrkas showed that (|7) viewed as a function of 
^, coincides the Chalykh-Veselov i/^-function (see [|], and Q for the matrix 
case considered in this paper), which is defined in terms of its behavior as a 
function of ^. 



3. The first formula 

Our first formula is an explicit calculation of the trace of the homomor- 
phism ^ : M^_p — > M^_p iSi U of the previous section. The image of the 
generating vector is a singular vector of weight — /> (a vector of weight — p 
killed by e^, a G A+) and all singular vectors of weight ^ — p correspond to 
some homomorphism. Let fi = e_a, , i = 1, . . . ,r and set fi = fi^ • ■ ■ fi„^ for 
a multiindex / = (ii. . . . , Let v^^p (8> ti + fi'^^-p uj he a singular 
vector in M^_p [/ of weight i — p. Our formula gives an eigenfunction 
in terms of the coefficients uj. Let us remark that, for generic ^, such a 
singular vector is uniquely determined by its first coefficient u € f/[0]: the 
coefficients uj can be given rather explicitly in terms of u and the inverse 
Shapovalov matrix (see ||6|). 

Theorem 3.1. Let ^ € f) * and let f^-p (!S)U + f^v^-p ®ul be a singular 
vector of weight — p in M^^p ® U . Then the function ip{X) = e'^'^^^^^\u + 
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with 




X 



where aj is the cardinality of the set of m € {j, . . . ,p — 1} such that a{m + 
1) < cr{m), belongs to E{S,), and for generic ^ all functions in -E(^) are of 
this form. 

Example. Let Q = sl2. If f/ is irreducible, U[0] is one-dimensional if U 
has odd dimension and is zero otherwise. Let U he a 2s + 1 dimensional 
irreducible representation. Our formula reduces to ^(A) = e'^'^^^^^\u + 
X]f=i(i^f {X = Xi). The components ui of the singular vector are 
easily computed, and we get the formula 



We will prove the more general quantum version of Theorem 3.1 in section 



Note that Theorem completes the proof of Theorem |2.1| : the coeffi- 
cients are meromorphic functions. They seem to have poles on {Xj3 = 1} for 
general (3 € Q; however, these poles cancel, since the differential equation is 
regular there: 

Lemma 3.2. Let e G C and suppose that ip is a meromorphic solution of 
the differential equation Hip = etp on i) , whose poles belong to the union of 
the hyperplanes Hp^m = {A|/9(A) = m}, £ Q, m G Z. Then ip is regular 
except possibly on the hyperplanes Ha,m, with a € A, m G Z. 

Proof: Let /3 G Q — A, m G Z, and choose a system of affine coordinates 
zi, Zr on 1} so that Zr = l3{X) — m. If tp has a pole of order p > on H^^rn, 
then, in the vicinity of a generic point of Hp^rn-, ip = Zr^f{zi, . . . , Zr~i) + . . . 
with non- vanishing regular /. The leading term of Hip as comes 

from the Laplacian, as the potential term is regular on Hp^m and is equal to 
p{p+l){a, a)zr^~'^f which is non-zero for p > 0. If Hip = eip, then it follows 
that ip has a pole of order p + 2, a contradiction. We have shown that ip is 
regular on generic points of Hp^^. Therefore, on any bounded open subset 
V oi \) , the product of ■0 by a suitable finite product of factors a(A) — 
a G A, / G Z is holomorphic on the complement of a set of codimension 2, 
and therefore everywhere on V by Hartogs' theorem. □ 



This concludes the proof of Theorem 2.1. 



4. The second formula 

Our second formula is an integral representation. It is obtained as the 
trigonometric limit of the integral representation of solutions of the Kni- 
zhnik-Zamolodchikov-Bernard (KZB) equations on elliptic curves with one 
marked point. 
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One ingredient in the integrand is the rational function of 2n variables 

W{T, y) = f]( -] , Tn+i := 1. 

Now suppose that U is an irreducible highest weight representation with 
non-trivial zero- weight space. The highest weight of U is then in i.e., of 
the form A = Ylj "^j^ji fo'^ some non-negative integers Uj. Set n = Uj. It 
is convenient to introduce the associated "color" function c : {1, . . . ,n} — > 
{!,..., r}, the unique non-decreasing function with |c^^({j})| = rij, j = 
1, . . . ,n. For each permutation a € Sn, let Wa^c be the rational function of 
n + r variables Ti, . . . ,Tn, Xi, . . . , Xr 

Wa,ciT,X) = W{T^(i),. . . ,r^(„),Xc(CT(i)), . . . (10) 

The other ingredient is the many- valued function of Ti, . . . , T^, 

n 

j<i i=i 
We will consider integrals of ^'^ j^{T)Wf^^ciT, X)dTi ■ ■ ■ dT^ over connected 
components 7 of {t G (0, 7^ tj, {i 7^ j)}. These "hypergeometric inte- 

grals" are defined as (meromorphic) analytic continuation in the exponents 
of Ti — Tj, Tj, Tj — 1 in ^>^^ from a region in which the integral converges 
absolutely. We say that a hypergeometric integral exists if this analytic 
continuation is finite at the given value of the exponents. 

Theorem 4.1. Suppose that U is an irreducible highest weight module with 
highest weight A = njaj and highest weight vector v\. Set n = Y^nj and 
let c : {1, . . . ,n} {1, . . . ,r} be the unique non- decreasing function such 
that c~^{j} has nj elements, for all j = l,...,r. Fix a generic complex 
number k and C G f) o,lso generic. Then, for each connected component 7 of 
{t G (0, / t.j, {i / j)}, the integral 

V'^(A) = e2-«(^) [ ^1^{T) W<x,c(T,X(A)) dT,--- dT^ • ■ ■ /,(.(„))z;a, 

•^T aG5„ (11) 

exists and defines a function in E{^). Moreover, for each generic k. and 4, 
all functions in -E(^) can be represented in this way. 

In the rest of this section we prove this theorem. The existence of the 
integral follows from |17], Theorem 10.7.12. Indeed, the coefficients in (11) 
are linear combinations of integrals considered in [17| in the context of the 
Knizhnik-Zamolodchikov equation with two points. 

Clearly ip.y is of the form e^'^*^ times a meromorphic function of Xi, . . . , X^ 
with poles on the divisors X/^ = 1, /? G Q. To prove that belongs to E{^) 
it is therefore sufficient, thanks to Lemma |3.2| , to show that ijjy is an eigen- 
f unction of H. 
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This follows from the results of 0, which we now recall. 

Let us define, for q = exp(27rir) in the open unit disk, 6{x) = vr"^ sin(7rx)nj°(l- 
2gJ cos(7rx) + g^i) (in the notation of 0, 0{x) = 6i{x) / 9[{id)) and v{x) = 
— -^ \nO{x). The function v is doubly periodic with periods 1 and r. Then 
for any k € C — {0}, the KZB equation is a partial diff'erential equation for 
a function n(A, r) on x {r € C [ Imr > 0}, with values in C/[0]: 

Attik— = An — i)(a(A))eae_aii. (12) 

As g — > 0, t'(x) — > TT^ sin~^(7r2;), and the differential operator on the right- 
hand side converges to —H. 

Proposition 4.2. Suppose that u is a meromorphic solution of the KZB 
equation ( [7^ such that q~"-/'^'^u is a meromorphic function of X and q, \q\ < 
1, and, as (? — > 0, 

n(A,T)=<7'^/2-(V^(A) + 0(g)) 
Then ip obeys Hip = etp, with e = (27r)^a. 

The proof consists of comparing the leading coefficients in the expansion of 
both sides of the KZB equation in powers of g at g = 0. 

A source of solutions of ([T^ ) with the property described in the Proposi- 
tion is 10]. In that paper, we gave integral formulas for solutions. 

A class of solutions with asymptotic behavior as above is constructed as 
follows (see Q for more details). Let U be an irreducible highest weight 
module with highest weight A € Q+, and let as above c : {l,...,n} — 
{1, . . . , r} the non-decreasing function associated to A. 

Solutions are labeled by an element ^ of and a connected component 7 
of {t G (0, l)"|tj 7^ tj, {i 7^ j)}. Let be a choice of branch over 7 of the 
many-valued function 

n 

i<3 i=i (13) 

and let w be the meromorphic function of ti, . . . , t^, yi, . . . , y„ 

and for each permutation a Sn, let w^^c be the meromorphic function of 
n + r variables ti, . . . ,tn, xi, . . . , Xr 

Wa,c{t, X) = W{t„(^i) , . . . , , , • • • , Xc(a{n)))- (14) 

Consider the differential n-form with values in U[0] depending on A S , 

= Wa,cit,ai{X), . . . ,ar(A)) dti--- dtn fc(a{l)) ■ ■ ■ fc(a{n))VA- 
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Then the integral 

%f,A(tV(A,T) 



7 



exists (as analytic continuation in the exponents from a region of absolute 
convergence) and is a solution of the KZB equation. As r — > ioo, 

6{x) = sin(7rx)/7r + 0{q), 

and the function (|^) behaves as 

in terms of the exponential variables Tj = exp(— 27ritj). Let Xj = exp(— 27riaj(A)). 
The components of the differential form uj behave as 

Wa,c{t, x)dti ...dtn = W„,c{T, X)dTi . . . dTn + 0{q). 

Therefore, the solutions corresponding to cycles in iJ„(C^(r), £(^)) have 
the properties of the Proposition, with a = (^,^), and if) is as stated in the 



claim of Theorem 4.1. We obtain in the limit an integral representation 
for eigenfunctions as in the Theorem but with an integration domain 70- 
of the form < aig{T„(^i^) < • • • < arg(To-(„)) < 2-k for some a € Sn- This 
integration domain may be deformed, so that the corresponding integral can 
be written as linear combination of integrals over the domains 7o-' defined by 
1 > Tg-z^i) > • • • > T^/(„) > 0. If we choose ^ = ia^o-, for some .^0 such that 
(^0) ctj) > for all j, and let a tend to infinity, then the integral over 70- is 
equal to the integral over 70- (for suitable choice of branch and orientation) 
plus terms that tend to zero. It follows that, for generic ^ and all a G Sn, 
the integral over 70- can be expressed as a linear combination of integrals 
over 7ct' and is therefore also an eigenfunction. 



This completes the proof of the first part of Theorem 4.1 



We must still prove that, in the generic case, all eigenfucntions in E{Q 



admit such an integral representation. In view of Theorem |2.1| , it is suf- 
ficient to show that for every u € C/[0], there exists a cycle 7 such that 
g-27n5(A)^^ — > li as ^ J = 1, . . . , r. In other words, one must show 
that all vectors in C/[0] are of the form 



n-l 



/ n^^T) ^ (H -) T • • • fci.in))VKdT, ■ . . dTn. 

for some cycle 7. But this follows from the results in [T^], Theorem 12.5.5. 



5. The third formula 

This is a formula of the Bethe ansatz type. 

Theorem 5.1. Suppose U is an irreducible highest weight module with high- 
est weight A = ^jrijOj and highest weight vector v/\_. Set n = ^rij and 
let c : {1, . . . ,n} {1, . . . ,r} be the unique non- decreasing function such 
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that c^^{j} has nj elements, for all j = I, ... , r. Then the function 
parametrized hy T ^ C" 

V'lr, A) = e2-«(^) ^(A))/c(a(l)) • • • fc{.{n))VK (15) 

(see ( [I^ for the definition of W^^c) belongs to -B(^) if the parameters Ti, 

. . . , Tn are a solution of the set of n algebraic equations ("Bethe ansatz 
equations ") 



Remark. These Bethe ansatz equations are the same as the Bethe ansatz 
equations of a special case of the Gaudin model (cf. 10). The solutions are 
the critical points of the function 

n 

{Ti - T,-)("<=W'"=(^)) J|rr^^"^'""<"'-'(rj - i)"(^'"'=o)). 

l<j<j<n j=l 



Theorem 5T can be understood intuitively as a consequence of Theorem 4.1 : 

one calculates the integrals, which up to normalization are independent of 
K, in the limit k — > 0, using the saddle point method. 

The proof of this theorem can be taken directly from ||^ in the case of a 
degenerate elliptic curve. 

The above result motivates the notion of Hermite-Bethe variety. Let 
A = ^jrijOj. Without loss of generality, we may take nj > 0, for all 
j = 1, . . . ,r (if this condition is not fulfilled, we may pass to a subalgebra). 
If c : {1, . . . , n} — > {1, . . . , r} is the associated non-decreasing function, we 
let Sc be the product of symmetric groups 3^ x • • • x Sn^- It acts on C" by 
permutation of the variables Tj with same c{j) and is a symmetry group of 
the system of Bethe ansatz equations. Moreover, we have, for all a £ Sc, 

V'(ar,A) =V(T,A). 

Let us write the Bethe ansatz equations as Bj{T) = 0, with 

Subtracting pairs of equations we may eliminate ^: 

Definition: Let F„ = (C - {0,1})" - yJi<r.^o.^^^^,a^^^^)^o{T\Ti = T,}. The 
Hermite-Bethe variety HB{c) is 

HB{c) = {T Fn\Bj{T) = B,+i{T),3 G {1, . . . , n}-{ni, m+na, . . . , n}}/5,. 

The remaining equations define a regular map p : HB{c) ^ , T i— > ^. 
The completeness hypothesis of Bethe states is in this case: 
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Conjecture 5.2. (cf. [18|, and conjectures in |T^^ The map p has 

dense image, and the generic fiber consists of (at least) dim(C/[0]) points. 
The eigenf unctions corresponding to points inp^^{^) span E{^). 



In the case of s/2 the conjecture holds and goes back to Hermite, see p9|| . 
We prove this conjecture in two special situations, including the scalar 
case relevant for many-body systems. 



Theorem 5.3. Conjecture |5.^ holds in the following cases: (a) q = sIn, 
U = SP^C^, p = l,2,..., (see (|) j (b) g = sIn, U = adjoint representa- 
tion. 

Proof: The proof is by induction in A^, viewing sIm as a Lie subalgebra 
of sIn+i- We choose simple roots ai,...,aAr of sIn in such a way that 
ai, . . . ,aN-i are simple roots of sIm-i- Also, it is convenient to replace ^ 
by C = ^ - 

In case (a), the highest weight of Un = S^^C^ is A = pYlf=i^ j'^j^ ^^'^ 
Um[0] is one-dimensional. The fiber over ^ = ^ -|- p of the Hermite-Bethe 
variety consists of the critical points of 

j<i j 

viewed as a function of T = (^i, . . . , Tp;Y(;Y+i)/2)- We have c(j) = m if 
pm{m — l)/2 < c(j) < pm(m + l)/2, 1 < m < A^. We prove inductively 
that, for generic C) ^^A^ has a non-degenerate critical point, and that the 
corresponding eigenfunction spans E{C + p)- Let us consider critical points 
of $Ar+i(C,T) when (ajv,C) = tends to infinity, and (ajX), j < ^ are 
kept fixed. Let ^ denote the orthogonal projection of ^ onto the Cartan 
subalgebra of sljy. Let us replace the coordinates Tj indexed by j such that 
c{j) = N hy new coordinates aj defined by Tj = 1 — eaj {c{j) = N), and let 
T = (Ti, . . . , Tpjv{Ar+i)/2) denote the remaining coordinates. Then 

c(j)=Ar-l c{j)=N~l,c{l)=N 
c{l)=N j<l 

This function converges as e ^ to a constant times 

^NiCf)lle^%'^''^'^ll{a,-a,f. 

The factor depending on the aj's has a non-degenerate critical point a^, in 
the domain Oj 7^ 7^ 0, (j 7^ k). This follows from the fact that this factor 
is obtained limit M ^ 00 of 

j 3<l 
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The critical points of the latter function are known explicitly (see (1.3.1) 
and (1.4.2) in [|l^]), and it can be easily checked that they have a limit as 
M ^ oo in the domain aj 7^ 7^ 0, which is therefore a critical point of the 
limiting function. In fact, a* is the set of zeros of the polynomial solutions 
of the differential equation xy" + (x — m)y' — Ny = 0, where m = p{N + 1). 
By the induction hypothesis, <I>Ar(C, T) has a non-degenerate critical point 
T!f. Moreover, in a neighborhood of (a*,T'*,e = 0), the right-hand side of 
(|l6|) is holomorphic, which implies that the non-degenerate critical point at 
e = deforms to a non-degenerate critical point for generic e. 

It remains to show that the corresponding eigenfunction generically spans 
the one-dimensional vector space S((^ -|- p), i.e., that it does not vanish. 
Assuming this inductively to hold for sljy, we see that as e ^ 0, the 
leading contribution in the sum ([l5|) is given by permutations such that 
c{a{j)) = N, whenever c(j) = N. These permutations give terms propor- 
tional to fc{T{i)) ■ ■ ■ fc{T{pN{N+i)/2)){fNT^VA, for some r G SpN{N+i)/2- It 
follows that when e ^ and T(e) is the above family of critical points, 
e-2'^*(«'^V^/iv+i(^(e), A) converges to C(a*)e-2'^*(?"^)^/;,ij^(f*, A), where the 
representation Un = S^^C^ of sIn is viewed as the s/Ar-sub module of Un+i 
generated by the singular vector {fN)^^VA- Using the identity 

\- I = ' 



we see that C{a^:) is, up to a trivial nonzero factor, the inverse of the product 
of the components aj of a*, and therefore nonzero. This completes the proof 
of part (a) of the theorem. 

The case (b) is treated in a similar way. The highest weight of the adjoint 
represnetation of sljsf is A = ^ ai and c = Id. The function $ for si jv is 

N-2 N-1 

i=i i=i 

As before we let e"-*^ = (Coat) go to infinity. If we set Tjy = 1 — ea, 
T = (Ti, . . . , Tjv-i), ( = projection of C onto the Cartan subalgebra of sIn, 
then 

I N~i — I + ea a 

Proceeding as before, we see that non-degenerate critical points T*, 0=,, = 1 of 
the limiting function deform to non-degenerate critical points for any generic 
e. the corresponding eigenfunctions deform to eigenfunctions of sIn taking 
values in C/Ar[0], viewed as a subspace of ?7Af+i[0] via the inclusion of sIn 
in s^Tv+i- In this way we get — 1 linearly independent eigenfunctions in 
the A^-dimensional space E{( + p) associated to the adjoint representation 
of sIn+1- To find the remaining eigenfunction, we let Tj = 1 — eaj for all j. 



14 



GIOVANNI FELDER AND ALEXANDER VARCHENKO 



Then as e — > 0, e'"^>Ar_|_i for suitable m converges to 

— — n " ^ 

The critical point of this function can be computed explicitly: aj = j(l + 
1/A^). It deforms to a non-degenerate critical point for e generic. The cor- 
responding eigenfunction converges, as e ^ 0, to a constant function. From 
its explicit expression it is clear that its value is not in Un = U{sl]\f)fNVA 
and is therefore linearly independent from the — 1 constructed before. □ 
Remark. The proof of the preceding theorem indicates, at least in the 
examples considered, that the construction which to a pair (g , U) consisting 
of a semisimple Lie algebra q and a finite dimensional q -module associates 
the closure of the algebraic Bethe-Hermite variety X is functorial: to each 
homomorphism {q,U) {g',U') preserving the Lie algebra and module 
structures there corresponds functorially a morphism of algebraic varieties 
X ^ X' . This functor is compatible with the construction of eigenfunctions 
in a sense that should be made more precise. 



6. WeYL group ACTION, JACOBI POLYNOMIALS, SCATTERING MATRICES 

Li this section we study the Weyl group action on eigenfunctions, and 
discuss the relation with multivariable Jacobi polynomials |8|. The Weyl 
groups W acts on U[0] (The normalizer of the Cartan torus T = expi) 
of the simply connected group with Lie algebra g acts on U, so W = N/T 
acts on f/[0]). Therefore we have a natural action of W on C/[0]-valued 
functions: w G W acts as {wip){X) = w ■'ip{w~^ X). The Schrodinger operator 
H commutes with this action. 

Let S be the set of G * such that (^, /3) = (/?, /?) for some (3 £ Q+. If 
^ G f) * — 5, E{^) is isomorphic to U[0] and we have the explicit expression 



of the isomorphism in Theorem 3.1 



Lemma 6.1. Let£,el)*-S. Ifip G and w€W, then e'^'"'^'^'" '^^u-V 
is a rational function of Xi, . . . ,Xr which is holomorphic on the complement 
in C of the root hypersurfaces X^ = 1, a G A. 

Proof: Suppose first that w = 1. From the explicit expression (|8|) it is clear 
that it is a rational function. By Lemma the poles lie on the root hy- 
persurfaces. Let us assume that u; = is a simple reflection corresponding 
to the simple root a^. As root hypersurfaces are permuted under the ac- 
tion of the Weyl group, it suffices to show that e"-^'^*^^"' ^^wtp is regular at 
Xj = for all j. This is obvious if j 7^ A; since replacing A by s^A amounts 
to replacing Xj by XjXj^ , with < if j ^ k. If j = k, X^. is replaced 
by XJ^^ and we have to check that the expression in parenthesis in (^) is 
regular when Xk — > 00, but this follows easily by counting powers of Xk in 
the numerator and denominator, exploiting the fact that aj + 1 < p — j + 1- 
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The case of general w is reduced to this case by writing w as a product of 
simple reflections. □ 

Lemma 6.2. Let G f) * 6e generic. For all w G W , ip i— > wij: is an 

isomorphism from -E(^) onto E{wS,). 

Proof: Let i/j € E{^). By the previous lemma, e"^'^*^^"' is holomor- 

phic at X = 0. Also, since H is invariant, wtp is an eigenfunction with 
eigenvalue 47r^(^,^) and therefore wil) G EiwS). □ 

Example. Let g =5/2, and U be the 2s + 1-dimensional irreducible represen- 
tation. The Weyl reflection si acts as (—1)* on C/[0]. Fix a nonzero u G U[Q\ 
and let be the eigenfunction @. Then (si'i/')(A) = (— 1)'''(/'^(— A) = 
5(^)'0-5(A), where S{^) can be computed in the limit X — > 0. This gives 
the "two particle scattering matrix" 

^ic^ ( 1^s (s + /)!r((g,a) - I) ^ fc + (a,0 

We say that a function -i/' is Weyl antiinvariant if w'i/' = e{w)tlj for all 
It; G T^. 

Proposition 6.3. (cf. Prop. 3 in ^) Let ^j) be a meromorphic Weyl anti- 
invariant solution of the eigenvalue problem Hip = eip, regular on 

[) - Umez UagA {A G I a(A) = m}, 

such that ip{X+p) = x(p)V'(A) for allp in the lattice = G f) | a{p) G Z} 
and some character x of P"^ , then tp extends to a holomorphic function on 
\] . Moreover, for all a ^ IS. and m ^ Z, 

eliP = 0{{a{X)-my+'), (18) 

as a(A) —f m. 

Proof: By periodicity with respect to the coweight lattice P^, we may limit 
our considerations to the hyperplanes through the origin. One proceeds 



as in Lemma 3.2 by approaching the singular hyperplane in a transversal 
direction. The leading term in the Laurent expansion in the transversal 
coordinate x = a{X) is determined by the differential equation 

{a,a)—Tipo 5 ipo = 0. 

Let us decompose U into irreducible representations of the subalgebra gen- 
erated by e±a. Since tp is of zero weight, we may replace CaS-a + e-aSa 
by ^ 

(a, a) 

where h^ = [ea,e^a\- If V'o belongs to a 2Z + 1 dimensional irreducible 
representation, i.e., if e^V 7^ but e^pp^ip = 0, then the Casimir element Ca 
acts as (a, a)l{l + 1). Therefore either ipQ ~ x'"*"^ or ipQ ~ x~K On the other 
hand, the Weyl reflection with respect to the hyperplane a = changes 
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the sign of x and multiplies the value by (—I)'- Since ip is antiinvariant, 
it follows that the first possibility is realized and the function vanishes to 
order / + 1. □ 

In particular, if ^ belongs to the weight lattice P, spanned by the fun- 
damental weights LOi, . . . ,ujr, and ip G E{^), then the Weyl antiinvariant 
function = X^^^iy is a rational function of X^^, . . . ,Xi^^ reg- 

ular on (C — {0})^. Therefore, V is a Laurent polynomial. To apply the 
previous lemma, ^ should not be in S. This is true if — ^ is dominant, and 
we have the following result. 



Corollary 6.4. Let ^ = —/j, where ^ is a dominant integral weight and tp G 
E{S). Then ip^ = '}2,we_w ^{w)w'ip is a Laurent polynomial in X^-^, . . . , X^,, . 

Example. Let g = sIn, U = S^^C^ , the scalar case, and let loi, . . . , wjv-i 
be the fundamental weights of sl^. Let us fix an identification of U[0] with 
C. Let V'-^ be the eigenfunction of the previous corollary, normalized so that 
^^g-27riMA) + By the previous corollary, it is a Laurent polynomial 



in the X^.. Then the vanishing property (18) of Proposition 6^ implies that 
the Weyl antiinvariant eigenfunction is divisible in the ring C[X^,^] by 
n^'+-^ where 

n = x„p J] (1 - x„) 

The ratio P = ip-^/W^^ is a Weyl invariant polynomial (since 11 is antin- 
variant) and obeys the differential equation 

-AP+{p + 1)2tt cot{7ra{X))daP = IP (19) 

aeA+ 

{da is the derivative in the direction of a), with e = A'iT'^{fi + {p + l)p, fj, — 
{p + l)/o)- This equation defines multivariable Jacobi polynomials (or Jack 
polynomials) associated to a dominant integral weight z/: The Jacobi poly- 
nomial Pi, is the unique solution in C[X^^, . . . , X^^__^] of ([l9|), normalized in 
such a way that P^ = X^^ + Yl,peQ+ ^p-^-u+P ^^^^ some coefficients aa G C. 
The leading term of ip-^/W^^"^ is X^_|_(p4.x)p- After antisymmetrization, we 
obtain a leading term cX^jj(-^4.(p+i)p) where wq is the longest element of the 
Weyl group. Thus we obtain a function proportional to P^ by choosing 
/i = wqv — {p + l)/5 {wq is an involution). Therefore we have the corollary: 

Corollary 6.5. In the scalar case g = sIn, U = S^^C^ , with fixed iden- 
tification of U[0] with C, let I' be a dominant integral weight, and for any 



^ G [) * — S*, let denote the eigenfunction in E{^) given in Theorem 3.1 
with u = 1 and ij^Y = X^^giy (-{w)w'ijj^ its antisymmetrization. Then 

p _ ^wou-{p+l)p 

for some constant Cj/ 7^ 0. Here wq denotes the longest element of the Weyl 
group, i.e., the permutation j N + 1 — j of Sn- 
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A formula for Cy is given below. 

We next study more closely the action of the Weyl group. 

Let (j) be the subalgebra of g generated by e±Qj. , and for a g -module 

f7, let U = ®sUs^^ be the decomposition of U viewed as s/20)-module into 
isotypic components: Us is isomorphic to a direct sum of 2s+l-dimensional 
irreducible sZ2(j)-niodules. 

Theorem 6.6. Let^ he generic and denote by the isomorphism U[0] — > 
E(^) mapping u to the eigenfunction with leading term e'^'^^^^^^u. Then there 
exists a family of maps S{^) : W GL{U[0]), w i— > such that if 

ip = then wil) = j{w^)Sw{C)u. These maps have the composition 

property 

If w = Sj is a simple reflection then Sw = Sj has the form 

where p]''^ G End(C/[0]) is the projection onto Ug''^ n U[0]. 

Proof: The first part of this theorem is just a rephrasing of the preceding 
lemma. The composition property follows from the definition: 

j{wiW20SwiW2{0 =WiW23{0 = Wij{w20Sw2{0 = j{wiW2i)Sw^{w2i)Sw^{i). 

If u) = Sj, we may compute Sj by letting all Xi, I ^ j go to zero. Then 
the coefficients Al in Theorem |3.1| tend to zero except if L = ... , j). 
and the formula for exp(— 27ri^(A))t/' reduces to an s/2 formula, and we can 
apply the previous example in each isotypic component. □ 
Remark. Note the similarity between this scattering matrix and the ratio- 



nal i?-matrix associated with general represetations of ,3/2 

Example. Let g = sIn, and identify f) * with C^/C(l, . . . , 1). The Weyl 
group is the symmetric group Sn and the simple reflections sj act on 
by transposition of the jth and j + 1st coordinates. The corresponding Sj 
depends only on (aj,^) = — (,j+i- The relations SjSj+iSj = Sj+iSjSj+i 
translate into "unitarity" 

Sji(j+i - Cj)Sj{^j - Cj+i) = Id, 
and the Yang-Baxter equation 

Sj {^j+i-Cj+2 )Sj+i {ij-ij+2 ) Sj i^j-<j+i) =Sj+i i^j-<j+i) Sj {^j-(,j+2)Sj+i{^j+i-Cj 

Corollary 6.7. In the scalar case (q = sIn, U = S^^C^ ) 



18 



GIOVANNI FELDER AND ALEXANDER VARCHENKO 



(7) 

Proof: The only isotypic component Us with non-trivial intersection with 
U[0] has s = p m this case. □ 



Let us conclude with a formula for the constant Cj/ in Corollary 6.5. This 
constant can be computed from the leading term in which appears in the 
summand indexed by wq: = e(wo)5'^o(^o('^ + + If we identify 

u £ i) with the diagonal traceless matrix with diagonal entries ui, . . . , i/^Vj a 
simple calculation gives 



-1 _ TT -rr k + ui-uj-{p + l){i-j) 



This constant is clearly different from zero since for dominant v we have 
Ui < Uj \i i > j. 

7. BeTHE ANSATZ IN THE ELLIPTIC CASE 

Let us fix r in the upper half-plane, and denote Er the torus C/Z + tZ. 
We first quote the result of M on eigenfunctions of the differential operator 



He = -A + ^ v{a{X))eae- 



on C/[0]-valued functions on f) . The elliptic function v is v{x) = —-^ ln9{x), 
6{x) = iT~^ sin(7rx)nf^(l - 2q^ cos(7r2;) -|- q'^^), q = e^'^*'^. It differs from the 
Weierstrass p- function by a constant. 

Theorem 7.1. (^) Suppose U is an irreducible highest weight module with 
highest weight A = highest weight vector v\. Set n = '^nj 

and let c : {l,...,n} — > {l,...,r} be the unique non- decreasing function 
such that c^^{j} has Uj elements, for all j = 1, . . . , r. Then the function 
parametrized by t £ C" 

llj{t, A) = e^''*^^^) ^ W^^cit, ai(A), . . . , ar(A))/c(a(l)) • • • fc{a{:n))VK 

aGS„ (22) 

(see dl^ ) for the definition of Wa^c) is an eigenfunction of if the param- 
eters ti, . . . , tn are a solution of the set of n equations ("Bethe ansatz 
equations") 

^e'(ti-ti) 9'{tA ^ ^ \ 

2^ nr; . ^ »c{i) - jrrT^ + 27rze, a^O) =0, j = l,...,n. 

The corresponding eigenvalue e is 

e = 47r2(e,0 -47rz^S(ti,...,t„,r), 
S{ti,. . . ,tm,T) = '^{ac{i),ac(j))lne{ti-tj) -'^{A,ac(i))lnB{ti). 

i<j i 
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Remark. Solutions of the Bethe ansatz equations are critical points of 

i<j i 

in the domain $e(i) 7^ 0. 

As in the trigonometric case, we define the Hermite-Bethe variety by 
eliminating the spectral parameter ^ from the Bethe ansatz equations. The 
resulting equations are the n — r equations 

A ^(i. - ^0 ^(t,+i - to r^^'^ [0{t,) 9(i,+i)/''"^(^)V (24) 



where j runs over the set of indices such that c(j) = c(j + 1). 



Lemma 7.2. The left-hand side of each of the equations ( \24 ) is a doubly 
periodic function of each of its arguments ti. 



This lemma follows easily from the formulas 

e'{t + i) 6'{t) e'it + T) e'{t) 



2Tri, 



e{t + i) 9{t)' e{t + T) 9{t) 

taking into account the zero weight condition A = ^ f^j'^j- 

Therefore we may view the equations as algebraic equations on E'^ = 
X • • • X E-r, a product of elliptic curves. Moreover, we have an action of 
Sc = X • • • X Srir permuting ti with the same c(z), that maps solutions 
to solutions, and does not change ■0. Let D{c) be the set of t G E"^ such 
that ti = tj for some « / j with (ac(j)) ctc(j)) 7^ 0, or = for some i with 
(A, ac(j)) 7^ 0. On this set the Bethe ansatz equations are singular. 

Definition: With the notation of the Theorem, assume that Uj > for all 
j. The Hermite-Bethe variety HB{c) is the subvariety of {E!^ — D{c))/Sc 
defined by the equations (|23) . 



The remaining r equations determine ^ as a function of a solution t of (24). 



They can be chosen to be the equations (|23D with j = ni,ni + n2, . . . ,n. 
We see from this formula that if tj is replaced by tj + n + mr, then ^ is 
shifted by —mOc^j)- It is easy to see that these replacements do not change 
the eigenfunction ip. Therefore we have a map : HB[c) — > f) */<5 mapping 
t to ^, and HB{c) parametrizes eigenfunctions tp such that ^{\ + uo) = 

One would like to prove that "all" eigenfunctions are obtained in this way. 
In the general case not much is known, however in the scalar case and the 
case of the adjoint representation we have the following result: 

Theorem 7.3. Let g = sIn, U = S'^^C^ or the adjoint representation. 
Then for each generic ^ G f] there are dim(C/[0]) solutions t G C" of the 
Bethe ansatz equations The corresponding eigenfunctions obey + 

oj) = e^'^^^^^^il){\), cj G , and are linearly independent. 
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The proof of this theorem is essentially the same as in the trigonometric 



case (Theorem 5.3): in the case of s/2 one has Hermite's result, and one 
proceeds by induction in A'^. The point is that one only has to use the 
asymptotic behavior of the theta functions involved, and this is the same as 
in the trigonometric case. 

8. The (/-deformed case 

We give here the (7-deformed version of our first formula. Our conventions 
for quantum groups are the following. Let q he a complex number different 
from 0, 1 or -1. We fix a logarithm of q, so that q"" is defined for all a G C. We 
normalize the invariant bilinear form on g in such a way that {aj,ai) E Z, 
for all j,l E {1, . . . , r}. The Drinfeld-Jimbo quantum universal enveloping 
algebra UqQ , a Hopf algebra, is the algebra with unit over C generated by 
fj, Cj, commuting elements kj, and their inverses kj^ , j = l,...,r with 
relations kjCikJ^ = g(°i'°')e,, kjfikj' = (?-("^'"')/,, ejfi - he, = 5ji{k, - 
^J^)/{q ~ 9~^)i deformed Serre relations Sa{q) = 0, a = 1, . . . , 2m, (see 
[Q]). The coproduct is defined on generators to be ^{fj) = /j 1 + k~^ fj, 
A{ej) = ej fcj + 1 (81 Cj and A{kf^) = kf^ kf^. We consider modules M 
over UqQ admitting a weight decomposition into finite dimensional weight 
spaces -/Vf[/i], fj, G f)*, on which kj acts as q(°':i'^^\ In particular, the Verma 
module of weight /i is the quotient = UgQ //(/u) by the left ideal /(//) 
generated by kj — ej, j = 1, . . . ,r. It is generated by its highest 

weight vector v^, the class of 1. 

If [/ is a finite dimensional UqQ -module, and <I> a homomorphism of UqQ - 
modules $ : M^-p ^£,-p ® U, we may define, as in the classical case, 

V e^'^*^('^)tr»^ r i<I) 

Recall that Xj = g-27riaj(A)^ -^^ classical case, we have for generic ^, 
q, and for each u in the zero- weight space U[0] of a finite dimensional UqQ- 
module U, a unique homomorphism of C/g0 -modules ^ : M^-p — > M^-p^U , 
sending the generating vector v^^p to v^^p u So, in the generic case, 
the "(/j-function is uniquely determined by g,^ and a vector u (z U[0]. It 
is a formal power series, but its explicit expression below shows that it is 
actually a meromorphic function of A. 

If g = sIn and [/ is a deformation of the pNth symmetric power of C^, it 
was shown by Etingof and Kirillov Q that ipW is a common eigenfunction 
of a commuting family of difference operators (related by conjugation by a 
known function to the ^Ar_i-Macdonald operators). 

As in the classial case, the image of the generating vector is a singular 
vector (a vector killed by Cj, i = 1, . . . , r) of weight — p and all singular vec- 
tors of weight C~ P correspond to some homomorphism. Set /l = /zj • • • 
for a multiindex L = [li. . . . , Im)- Our formula gives, for each u E C^[0], "0 in 
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terms of the coefficients n^, in the formula of the the unique singular vector 
of the form v <Si u + fhv^-p <^ ul These coefficients are given explicitly 
in terms of u and the inverse Shapovalov matrix, see |p. 

Theorem 8.1. Lei ^ G * be generic and letv^-p^u+Y^^^ Jlv^-p^ul be a 
singular vector of weight p in M^_p ® U. Let <I> G Homfy^g {M^_p, M^^p ® 
U) be the corresponding homomorphism, and ip the ip-function (|^^. Then 
^{\) = e2^^«W(u + ^^^L(AK), with 

(p -^aj{a)+l \ 

where aj{cr) is the cardinality of the set of m G {j, . . . ,p — 1} such that 
a{m + 1) < cr{m), Xj = exp(— 27riaj(A)), and 

p mm 
J=l k<j,a{k)>a{j) meS j=l j=l 

S = {m e {1, . . . ,p - l}\a{m) > a{m + I)} U {p}. 

We conclude this section by giving a conjectural formula for Macdonald 
polynomials, which is a g-deformation of Corollary The ^at-i Mac- 
donald polynomials |12] are symmetric polynomials Pu{x,q,t) (we use the 
notation of Q) in variables xi, ... xjm, depending on parameters q and 
t. They are labeled by dominant integral weights i' of gljsr, i.e., decreas- 
ing sequence ui > ■ ■ ■ > of non-negative integers. We consider the 
case where t = for some positive integer A; as in Q. The symmetric 
polynomials P{x, q, q^) are uniquely characterized by having leading term 
■ ■ ■ x'^ (as Xi/xi+i — > oo, i = 1, . . . , — 1) with unit coefficient and by 
being orthogonal with respect to the inner product {f,g) = constant term 
of f{xi, . . .,XN)gix^^, . . . x^^)A(x), where 

k 

A{x)=YlYl{l-q'^x^/x,). 

i-^j m=0 

Obviously -Py+(i,...,i) = x\. . . xjyPi, so it is sufficient to evaluate Macdonald 
polynomials on the hypersurface xi ■ ■ ■ x^ = 1- Then v may be condidered 
modulo Z(l, . . . , 1), i.e., as sIn dominant weight. The (^-analogue of Corol- 
lary |6.5| is the following conjecture, which can be proved with the same 
method as in the classical case if = 2 or 3, but is open in the general case. 

Conjecture 8.2. Consider the scalar case g = sl^, with = {A G | ^ Aj 
0}, U = SP^C^, and fix an identification of U[0] with C. Set Xj = e^"^^ . 
Let 1/ G f) * be a dominant integral weight, and for ^ ^ , let ip^ denote the 



function given in Theorem p. j| with u = 1 Then 



n{x,q,q j -c, 2^ 
w&Sn 
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where 

p 

m=0 j<l 

for some constant Cu ^ 0. Here wq denotes the longest element of the Weyl 
group, i.e., the permutation j N + 1 — j of S^. 

As before, the constant may be computed in terms of the (^-analogue 
of the scattering matrix ( ^T|) 

x\ ^ 



^ [m + i^i- Vj -{p + l){i - j)]q _ q"" -q 

i>j m=0 



UU [rn-u, + u, + {p+l){i-j)],' '-^"^ 



9. Proof of Theorems 3.1 and 8.1 



We give the proof of Theorem |8.l| . The proof of Theorem ^Jj is essentially 
a special case. 

The proof is in two parts. We first prove that the calculation can be 
done in the algebra without Serre relation (Prop. |9.2D , and then do this 
calculation, reducing it to a combinatorial problem. 

We denote by Uqh the subalgebra of UqQ generated by fj, k^^, j = 1, . . . ,r 
and by UgXi the subalgebra generated hy fj, j = 1, ... , r. We have A(C/qn) C 
Uqb ^ UgXi. Let us introduce a Q-grading of these algebras by setting 
deg(/j) = aj, deg{kj) = 0. What we need to know about the deformed 
Serre relations is that (i) UgXi is the quotient of the free algebra generated 
by /i, . . . , /r by the ideal Jg generated by the deformed Serre relations si(g), 
. . . ,Smiq), (ii) si{q), . . . Sm{q) are homogeneous polynomials in the fj with 
coefhcients in Z[q,q~^] reducing to the (classical) Serre relations at g = 1, 
and (iii) A{sa{q)) = Sa{q) ^ 1 + K (g) Sa{q) for some K = Y[kj"^\ 

We will also need the fact that the dimensions of the spaces of fixed degree 
in Uqn are independent of q and coincide with the dimensions of the classical 
enveloping algebra C/n of the Lie subalgebra n C g generated by fj = c-aj , 
j = l,...,r (see, e.g., 

The first observation is that the computation of the trace, given the com- 
ponents of the singular vector is a computation in UqXi: 

Lemma 9.1. Let for w G M^^p, <I>"(u;) : M^_p M^^p^UqU he the unique 
homomorphism of UqXi-modules mapping v^^p to w<Sil. Then, in the notation 



of Theorem \8A 



^ E,e2-MA)tr,,_^[,]l 

This lemma follows immediately from the fact that Verma modules are free 
over UqXi. 

Let us introduce algebras without Serre relations. 
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Definition: Let Uqb be the algebra over C generated by /i, . . . , /r, kf^, . . . , k^^ 
with relations kjf^kj^ = q~^"^'°'''^ fk, kjki = kikj, kjkj^ = 1 and with 
coproduct as in UgQ . Let Ugfi C Ugb he the free algebra on generators 
fi, . . . , fr. The Verma module over Uq b of highest weight /x € * is the 
left module Uqb/I{pL) where /(//) is the left ideal generated by kj — 
J = 1, . . . , r. The image of 1 in is denoted by v^. 

The Verma module has a weight decomposition with finite dimensional 
weight spaces, and is freely generated by [/gfi-module. In particular, 

= M^/JqV^. Moreover, A(C/qn) C Uqb'^Uqn. Therefore the construction 
of the homomorphisms ^^{w) : M^-p — > M^-p ® Uqn works also in this case 
and traces make sense (as a formal power series). 

Proposition 9.2. Let for w G M^„p, ^"{w) : M^_p M^_p (g) UqX\ he the 

unique homomorphism of UqXi-modules mapping v^^p to w (8> 1, and for w G 
M^_p let ^^{w) be the same object for the algebra without Serre relations. 
Then for any w projecting to w under the canonical projection M^^p — s- 
M^-p, 



(26) 



The proof of this proposition requires some preparation, and will be com- 



pleted after Lemma 9.7 below. 

Lemma 9.3. If x € Uqn is homogeneous of degree a = rrijaj € Q 

Ij kj"^' (g)x + ^ Xj (g) Xj, wheri. 



then A{x) can be written as A(x) = J| • A; ■ '"^ (g x + ^ x' x", where x'' are 



homogeneous of degree < a. 

This lemma is an easy consequence of the form of the coproduct of the 
generators fj. 

Definition: Let n be the free Lie algebra on r generators /i,... ,/r (see 
[]lO| ). An iterated bracket of length one is an element fj of this set. An 
iterated bracket of length / > 1 is defined recursively as an expression of the 
form [a, b], where a, and b are iterated brackets of length li, l — li respectively, 
for some /iG {!,...,/ — 1}. An iterated bracket of length I is called simple 
if it is of the form [fi-,,[fi^, [• • • , /jj]]. 

Lemma 9.4. Let a,b G fi, and a be an iterated bracket. Then [a, b] is a 
linear combination of elements of the form 

[fn,[f..,A--- AkM]- (27) 

Proof: We use induction in the length of a. If the length of a is one there 
is nothing to prove. If a = [01,02] is of length I, with of length li < I, 
the Jacobi identity gives [a,b] = [ai, [02,6]] — [02, [ai,b]]. By the induction 
hypothesis, b' = [ai,b] and b" = [02, b] are a linear combination of elements 
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of the desired form (27). Then we use the induction hypothesis once more 
with b replaced by b' and b" . □ 

We apply this lemma to the following situation. The Lie algebra n C 
is the quotient of the free Lie algebra fi on r generators by the Lie ideal s 
generated by Serre relations si, ... , Sm, which are simple iterated brackets 
in it. 

Lemma 9.5. Denote by fi the free Lie algebra on r generators fi, ■ ■ ■ ,fr, 
and let si,...,Sm be simple iterated brackets in n. Let s be the Lie ideal 
generated by si,. . . ,Sm, and assume that d = dim(n/s) < oo. Then there 
exists a basis 61,62, •• • of n consisting of simple iterated brackets such that 
t>d+i^ i>d+2-, ■ ■ ■ is a basis of 5, and such that, for j > d, bj is of the form 

[/iiJ/iaJ--- J/ipSfc]]]]- 

Proof: By definition every element in fi is a linear combination of iterated 
brackets. The ideal s is spanned by iterated brackets containing at least 
one of the Sj. Using the skew-symmetry of the bracket, we see that every 
element of s is a linear combination of elements of the form [ai, [02, [. • • , Sfc]]] 
for some iterated brackets aj. The claim then follows from the preceding 
lemma. □ 

By the Poincare-Birkhoff-Witt theorem, the universal enveloping algebra 
Un has a basis b^ = IP^ ■ ■ - hi^y^ labeled by the set J of sequences J = 
(ji)i2i • • • ) of non- negative integers with ji = for all sufficiently large /. 
As the Verma module is freely generated by the highest weight vector 
v^j_ as a f7ti- module, {h''v^)j^j is a basis of M^. 

Let us extend this to the g-deformed case. Note that UqXi as an algebra 
is the same for all q. 

Lemma 9.6. There exists a sequence 61(g), 62 ((z), •• • of homogeneous ele- 
ments in Un ® C[q, q^^\ such that 

(i) 6,-(g) =6, mod {q-l)C[q,q-^] 

(ii) If q is a transcendental number, the elements b'^{q) = ■ ■ ■ b2{qY^bi{qy^ 
form a basis of C/gfi. 

(iii) For j >d + l, A(6j(g)) = bj{q) 1 mod Uqb (g) Jq. 

Proof: We may take bj{q) = bj if j < d. If a G UqXi is a homogeneous 
element of degree /3 G Q, set adq{fj)a = fjO — q^^'"^^afj. li j > d and 
bj = ad(/jj • • • ad(/jjsfc, set 

bj{q) = adq{fj^) ■ ■ ■ adq{fjj)sk{q). 

It is clear that (i) holds, (ii) holds too, since the determinant of the matrix 
expressing, in each homogeneous component of Uh^C[q, q~^], b^ [q) in terms 
of 6^ has a determinant in Q[(7, g"^] with the value 1 at g = 1, and is 
therefore invertible if q is transcendental. As for (iii), we know that Sk{q) 
has the required property since A(s/;(g)) = Sk{q) I + K (g Sk{q) for some 
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K. Moreover if a € UgVi of degree /3 obeys A(a) = a® 1 mod Ugb ^ Jq, then 

A(adg(/j» = fja 1 + /c~^a fj - q^^'^^'Hafj (g) 1 + akj^ (S> fj) mod Uqb (g> Jq 

= adg(/j)a (8) 1 mod Uqb (^Jq, 

since ^j'^a = q^^'^'^^akj^ . Therefore, by induction, Abj{q) = bj{q) 1 
mod f7g6 ^ Jq. □ 



Suppose that g is transcendental, and fix a sequence bj{q), as in Lemma 9.6 
and thus a basis {b"^ (q)) of ?7qn. To simplify the notation, we will write, for 
any J £ J, deg(J) = deg(6'^(g)). Let J", be the set of J = (ji, j2, ■■■) & J 
such that = if / > d, and J"' be the set of J such that jz = if / < d. 
Then J = J' J" canonically, and we may write the basis as b"^ {q)b'^ (q), 
{J', J") e J' X J". Then bJ'{q)bJ"{q) G Jq if J' is non-trivial. Since the 
dimensions of the homogeneous components of UqW are the same as in the 
classical case, the basis elements with J' non-trivial form a basis of Jq and 
the classes of b"^" (q) form a basis of UqXi = Uqn/Jq. 

Lemma 9.7. Suppose that q G C is transcendental. Then the elements 
b-^^q), J ^ J' form a basis of the subalgebra of UqU consisting of elements x 
such that A(x) = x ®1 mod Uqb ® Jq 



Proof: Call B this subalgebra. By Lemma |9.6| , (iii), the linearly independent 
elements b'^{q), J ^ J' belong to B. What is left to prove is that any 
X & B can be written as a linear combination of these elements. Write x = 
J2j'xj" C'J',.J"^'^ {q)b'^ {q)j with complex coefficients aj'^j/'. Let Z C Q+ be 
the set of /? such that there exist J', J" with deg(J") = (3 and aji^jn ^ 0. 
An element /3 € Z is called maximal if /3' > /5 implies f3' Z. Obviously, 
every element in Z is < some maximal element. Therefore, if we show that 
the only maximal element is 0, we have proved the lemma. 
By Lemma I 



A(2;) = ^ aji^jn {b-^' (q) 1) A(6"'" (q)) mod Uqb (g Jq. 

Let P be maximal. The terms whose second factor has degree P are, by 
Lemma 

A(x) = h ^ aj>j"b-^'{q) (g) b-^"{q) H mod Uqb (g) Jq. 

dcg{J")=/3 

These terms must vanish if /3 ^ since x B. But the elements 6"^" are 
linearly independent, even mod Jq. And, by construction, aj'j'i ^ for 
some J', J" with deg( J") = /?. Therefore /3 = 0. □ 

The proof of Proposition |9.2| is a calculation of the traces using the basis 
b-^{q)v^-p, J e J of M^-p and b-^{q)v^-p, J G J" of M^^p. It is sufficient to 
prove Proposition free in the case where q is transcendental, since all traces 
over weight spaces are clearly rational functions of q. Let w and w be as 
in Proposition [9^ . If J G J", Ab-^{q)w = ^L^jb^{q) (g) Ajl, for some 
AjL G Uqh. 
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The numerator of the left-hand side of (1261) is then 



^e2-''WtrA,,_^[^]Cl>»= ^ AjjX^_^+deg(j) mod J,. 
^J. Jej" 

On the other hand, the calculation of the numerator on the right-hand side 
of (|26|) involves (see Lemma |9.6| , (iii)) 



A{b-^' {q)b-^" {q))w (S)1=Y^ h-^'{q)h^{q) ® Aj.l mod Ugb Jg. 

Lej 

Now we claim that the only terms in this summation that may give a non- 
trivial contribution to the trace are those with L G J'". Let indeed b^{q) = 
b^'{q)b^"{q), with L' G J' and L" G J". Since both b-^' (q) and b^' {q) are in 
the algebra defined in Lemma |9.7| , their product is also in this algebra and is 
thus a hnear combination of b^q), M' G J', deg(M') = deg( J') + deg(L')- 
If L' ^ 0, this degree is strictly larger than deg(J') and therefore does not 
contribute to the trace. The trace is therefore 

J^e'"''^^^tr^,_,M^"(*) = ^ Yl ^J"J"^P-€+dcg(J")^deg{j') mod J, 
Ai J'eJ' J"ej" 

= (E ^deg(J'))E^''^^'^'^t^M.-pM^"H 
J'&J' 

On the other hand, 

^g27rt/.(A)^j,^^^^^^^-^ = E E -^P+deg(J")+dog(J') 

M J'ej' J"ej" 

Taking the ratio of the last two expression completes the proof of Proposition 



9.2 



We now turn to the actual calculation of the ratio of traces on the right- 
hand side of (|2^). We use the basis (/j = /j^ . . . fj^) of Uqn, labeled by 
finite sequences J in {1, ■ ■ ■ ,r} of arbitrary length m > (if m = 0, we set 

^ = 1), and we may assume that Hj = fivs^-p, with L = (/i, . . . , Ip). By 
definition, the trace in the numerator of (26) is the sum over J = (ji, . . . ,jm) 
of the diagonal elements Bjj in 

A{fj)fLVi-p ® 1 = E /a-^^^-P ® ^A'^^' (28) 
M 

weighted by a factor ^p-g IlfcLi -''^jfc • Expanding the coproduct on the left- 
hand side of (|2^) gives 2™ terms, each of them of the form 



(ffiH ■ ■ ■ fhsfLV^-p ® fu^ ■ ■ ■ fj, 
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labeled by subsets B = {bi < ••• < bs} of {l,...,m} with complement 
D = {di < • • • < dm-s}- The exponent of q is 

a = {C- P + ^ai^,^aj^)+ J] (ajM^id)- (29) 

j deD B3b>deD 

This term gives a contribution to the trace if and only if 

{jbi,---,jbs,h,---,lp) = ijl,---,jm)- (30) 

If it does, the contribution is 

q^X^^pXj^ ■ ■ ■ Xj^fj^^ ■ ■ ■ . 

Note that the j^- are necessarily equal to the k up to ordering. In par- 
ticular m — s = p. The terms contributing to the trace in the numera- 
tor are therefore in one-to-one correspondence with pairs consisting of a 
finite sequence (ji, . . . ,jm) and a subset C {1, . . . , m} obeying (|30[). Let 
S = {1 < 6i < 62 < • • • < 6s} C N be given. If bj = j for all j, set t = s + 1. 
Otherwise let t be the smallest integer so that bt > t. Then the condi- 
tion (|3^) on J determines uniquely jt,jt+i, ■ ■ ■ ,jm in terms of L, and gives 
no constraint on j'l, . . . ,jt-i- Therefore we can restrict our attention to the 
case t = 1: the general solution (J, B) of ( pO| ) is J = (ji, . . . ,jt-i,j'i, ■ ■ ■ ,j'm)y 
B = {1, . . . ,t - l,b[ + 1 - 1, . . . ,b'^ + t - 1} where ( J', B') is a solution with 
b'l > 1, and ji, . . . ,jt~i are arbitrary. The contribution to the trace of such 
a solution ( J, B) is • • • Xj^_-^ times the contribution of ( J', S'). The sum 
over all such solutions with fixed ( J', B') gives a factor that cancels with the 
trace in the denominator (except for Xp). 

00 

m=Oji,...,jm 

This implies 

Lemma 9.8. Lei L = {li, . . . ,lp). 



V p27rj^(A)tr - 1 ~ 'i ^Jl ^JmJ Jd^ ■ ■ ■ J 3d„^ 



where {J,B) runs over all solutions of ( p^ ) suc/i t/ioi 61 > 1. T/ie exponent 
a = a{J,B) is (^^, and {di < ■ ■ ■ < dp} is the complement of B. 



Our problem is therefore reduced to the combinatorial problem of finding all 
solutions {J,B) of ( |30| ) with 5i > 1 and computing their contribution to the 
trace. The problem can be reformulated as follows: consider a circle with 
p distinct marked points numbered counterclockwise from 1 to p starting 
from a base point * distinct from the marked points. The marked points 
are assigned labels: the jth point is assigned the label Ij. Solutions {J,B) 
are in one-to-one correspondence with games that a player can play erasing 
marked point on this circle according to the following rules. The player walks 
around the circle in clockwise direction starting from * a finite number of 
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times to return to *. Each move consists of proceeding to the next (not yet 
erased) marked point. When the player meets a marked point, he or she 
has the option of erasing it. If he or she does not erase it, the score is Xi 
where / is the label of the point, and we say that the player "visited" the 
point. The score for erasing the point is Xiq^°^^'^\ where f) = nijai- and 
nij is the number of times the jth. point has been visited previously. The 
game continues until all marked points have been erased. The score of the 
game is the product of scores collected during the game. For instance, if 
L = (1,1,2), one game consists of the sequence VEVEE of visits (V) and 
erasures (E). Its score is 

The relation with the previous description is that jm, ■ ■ ■ ,j2,ji is the list of 
labels of marked points met (visited or erased) at each move of the game, 
and B indicates at which moves points are visited. The contribution of a 
solution {J,B) to ( |3l| ) is the overall factor 

times the score of the corresponding game times fi^ - ■ ■ fii where ii, . . . ,ip 
are the labels of the erased points, in order of erasure. 

The next step is a reduction to the classification of "minimal" games. An 
empty round in a game is a sequence of of subsequent visits of all marked 
points exactly once without any erasures. A game without empty rounds is 
called minimal. Out of a game we can construct a new game by inserting an 
empty round just before an erasure. The score of the new game is HseE -^is 
times the score of of the old game, where the product is taken over the set E 
of marked points yet to be erased, and b = \ "^seE '^hl'^ ■ game can be 
obtained from a unique minimal game by doing this construction sufficiently 
many times. The games obtained this way have all the same sequence of 
erased points, and give thus proportional contributions to the trace. 

Minimal games are in one-to-one correspondence with permutation a G 
Sp-. we erase first a{p), then a{p — 1), and so on, always at the earliest 
opportunity. The example above gives the minimal game corresponding to 
the permutation (231). If a is the identity, the score of the corresponding 
minimal game is Xi-^ ■ ■ ■ Xi^. For general cr, we must calculate the numbers 
bjk of times the point a{k) has been visited before a{j) is erased. The score 
of the minimal game corresponding to a is then 

n " = ^^■'^("'.0) ' (33) 

3=1 j,k=l 

The sum of the scores of all games obtained from this game by inserting 
empty rounds is then 



(34) 
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We proceed to compute the numbers bjk- Roughly speaking, bjk is the 
number of times one goes around the circle before erasing either a{j) or 
a{k). More precisely, let 

5^ = {m e {!,... ,p- l}\cr{m) > a{m + 1)}. 

li k < j and cr^k) > cr{j), a{j) is erased first, and bj^ is the number of times 
the player passes through the base point *, including at the beginning of 
the game, before erasing (T{j). This number is bjk = \{m G S^\m > j} \ + 1. 
If A; < J, and a{k) < cr(j), bjk is the number of passages through * before 
erasing minus 1: bjk = \{m £ S^\m > j}\. li k > j, we get in all cases 
the number of passages through * before erasing cr{k), minus 1: bjk = \{m E 
S^\m > k}\. In particular, bjj = \{m G {j, . . . ,p — l}\a{m) > a{m + 1)}. 
The exponent c of g in (^) is then 

m m 

fc<j,o-(fc)>(T(j) me5x j=l i=l 

The formula in the theorem is then obtained by combining this formula with 

(H) , (H). 
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